Anti de Sitter black holes and branes in dynamical Chern-Simons gravity: 
perturbations, stability and the hydrodynamic modes 
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Dynamical Chern-Simons (DCS) theory is an extension of General Relativity in which the grav- 
itational field is coupled to a scalar field through a parity violating term. We study perturbations 
of anti-de Sitter black holes and branes in such a theory, and show that the relevant equations 
reduce to a set of coupled ODEs which can be solved efficiently through a series expansion. We 
prove numerically that black holes and branes in DCS gravity are stable against gravitational and 
scalar perturbations in the entire parameter space. Furthermore, by applying the AdS/CFT dual- 
ity, we relate black hole perturbations to hydrodynamic quantities in the dual field theory, which 
is a (2-|-l)-dimensional isotropic fluid with broken spatial parity. The Chern-Simons term does not 
affect the entropy to viscosity ratio and the relaxation time, but instead quantities that enter the 
shear mode at order in the small momentum limit, for example the Hall viscosity and other 
quantities related to second and third order hydrodynamics. We provide explicit corrections to the 
gravitational hydrodynamic mode to first relevant order in the couplings. 



I. INTRODUCTION 

General Relativity (GR) stands as one of the most 
elegant physical theories, and it describes gravitational 
interactions from small to at least solar system scales 
with high accuracy. Nevertheless, there are several in- 
dications that Einstein's theory should be modified, this 
"evidence" ranging from unification to quantization argu- 
ments. A variety of different proposals to extend GR have 
been made ever since its conceptual formulation almost a 
century ago. A promising extension of GR is Dynamical 
Chern-Simons (DCS) gravity ^ , in which the Einstein- 
Hilbert action is modified by adding a parity-violating 
Chern-Simons (CS) term, which couples to gravity via 
a scalar field. This correction arises in many contexts. 
Such a term could help explaining several problems of 
cosmology, from inflation (as discussed by Weinberg Q) 
to baryon asymmetry . In most of the moduli space 
of string theory, a CS correction is required to preserve 
unitarity; furthermore, duality symmetries induce a CS 
term in all string theories with a Ramond-Ramond scalar 
0. In loop quantum gravity, the CS term is required to 
ensure gauge invariance of the Ashtekar variables [loj and 
it also arises naturally if the Barbero-Immirzi parameter 
is promoted to a field [ill 113 ■ 

DCS gravity is described by standard GR supple- 
mented with a parity violating term, which can be 
thought of as the parity violating counterpart to the 
Gauss-Bonnet term: in analogy to the Gauss Bonnet 
term, in four dimensions the CS term is a topological 
term, which enters the equations of motion only when 
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coupled to a scalar field. Due to the parity violating na- 
ture of the CS term, spherically- and plane-symmetric so- 
lutions in GR, such as the Schwarzschild or the Friedman- 
Robertson- Walker-Lemaitre spacetime, are still solutions 
of the DCS theory. Nevertheless, these spacetimes, which 
are simultaneously solutions in different gravity theories, 
respond differently in both theories to generic external 
perturbations: the perturbed spacetime needs not, and 
in general does not, obey the underlying symmetry of 
the background solution and therefore probes more com- 
pletely the structure of the field equations. This point, 
which was raised by many authors (see e.g. Ref. [l3j). 
can be clearly seen by considering perturbations of a 
Schwarzschild black hole in DCS theory: in GR, the 
linearized equations can be reduced to two decoupled 
second-order ODEs and both gravitational polarizations 
have the same characteristic spectra [31; in DCS grav- 
ity however, the axial sector of the theory couples to the 
scalar field, resulting in two coupled ODEs for this sector, 
and a loss of isospectrality in the oscillations, or quasi- 
normal modes (QNMs) [TallBl- In flat spacetime, non- 
rotating black holes in both theories are stable |la |. 

For asymptotically anti-de Sitter (AdS) spacetimes, 
most of the previous considerations still hold. However, 
so far investigations of DCS gravity have been almost 
only limited to the case of asymptotically flat spacetime, 
with the exception of Ref. [13] , where some solutions to 
the non-dynamical version of CS gravity with a nega- 
tive cosmological constant have been studied. Very re- 
cently, DCS gravity in asymptotically AdS spacetime has 
been considered in the context of gauge/gravity corre- 
spondence [l8i . However no study has yet been made of 
either the stability properties or the oscillation modes of 
AdS black holcs/branes in DCS gravity. 

AdS black hole solutions are of particular interest in 
relation to the celebrated AdS/CFT correspondence fl9| . 
This duality has been originally established between 
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d = 4, = 4 Yang-Mills theory and a string theory in 
^^5*5 X . However, it is often seen as a particular real- 
ization of a more general gauge/gravity duality, mapping 
AdS solutions of the gravity theory in d + 1 dimensions 
to conformal field theories living on the d dimensional 
boundary spacetime and viceversa. In this context, AdS 
black holes arc dual to thermal quantum field theories at 
strong coupling and their characteristic oscillations dom- 
inates the near-equilibrium response of the dual gauge 
theory. More precisely, black hole QNMs in the bulk 
are related to the poles of the retarded correlators of the 
corresponding dual operator, providing useful insights on 
the transport coefficients and on the quasiparticle spec- 
trum (see Refs. [13, H^l for a review). Within Einstein's 
gravity in five dimensions, black hole perturbations have 
been extensively investigated and have revealed unex- 
pected properties of the corresponding strongly coupled 
field theory. For example, the celebrated universality of 
the shear viscosity to entropy ratio [U [HI can be under- 
stood in terms of the universality of black brane QNMs 
in the hydrodynamic limit. 

Considering the holographic duals to alternative theo- 
ries of gravity is interesting for several reasons. First of 
all, any consistent correction (as the DCS term) to GR 
can in principle modify the dynamics of the dual the- 
ory, with potentially interesting effects. There has been 
an intense effort to apply this approach in order to de- 
scribe holographic duals of strongly coupled system from 
a phcnomenological point of view. Furthermore, adding 
second-order in curvature corrections to the gravity the- 
ory corresponds, from the holographic point of view, to 
consider the next to leading order corrections in a strong 
coupling expansion. 

Since Schwarzschild black holes are still solutions of 
DCS gravity, the holographic theory will be the same as 
in GR, but the CS coupling will affect the linear response 
of the dual theory at equilibrium. For example hydrody- 
namic quantities and quasi-particlc spectra will carry the 
imprint of the CS parity violation. Very recently [isj . 
DCS gravity has been considered as a simple holographic 
realization of a (2 -I- l)-dimensional isotropic fluid with 
broken spatial parity in a bottom-up approach, and it 
was shown that the holographic theory may possess a 
non-zero Hall viscosity. The Hall viscosity is related to 
the presence of a non-trivial background scalar field, en- 
coding the nature of the parity violation in CS gravity. 

In this work we shall follow the same bottom-up 
approach and investigate the perturbations spectrum 
of Schwarzschild-anti-dc Sitter (SAdS) black holes and 
branes in four dimensions. The (2 + l)-dimensional dual 
field theory would be defined on a two dimensional sphere 
(or plane, for black branes) and it will be free of trian- 
gle anomalies [22|, ■ Because the SAdS back grou nd is 
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van- 



even-parity, the Hall viscosity studied in Ref. 
ishes, at least at first order. Nevertheless, the CS cou- 
pling will affect the hydrodynamic limit of the dual the- 
ory, in a characteristic way which we explicitly compute. 
Our main results are that 



(i) black holes and branes in this theory are stable. 

(ii) in the hydrodynamical limit, the shear mode acquires 
DCS corrections. We find 



^shcar 



where 7 = ^-^^osS+x/stt _ q_q28i and the novel correc- 
tion factor X is 



X= g^(201 - 20\/3^-601og3 



0.123072. 



(2) 



Here, a is a DCS coupling (defined below), r^, L arc the 
horizon and the AdS radius respectively, 
(iii) the sound mode is unaffected. In particular, the 
ratio rj/s = l/47r still holds in DCS gravity. 



II. SETUP 

DCS gravity with a negative cosmological constant is 
described by the following action 




d^Xy^Va(f'^''(, 



(3) 



where (/) is a real scalar field, which couples to the Pon 
tryagin density 

1 
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Without loss of generality, we set k = 1 = /3 and write the 
cosmological constant as A = — where L is the AdS 
curvature radius. The variation of the action with 
respect to the scalar field and to the metric component 
leads to the following system of coupled equations: 



Rab — —aCab 
a 
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RR. 



where 



(5) 
(6) 

(7) 



We consider spherically- and plane-symmetric black 
objects, which we write in a unified form as [23 | 

ds2 = g^^dx'^dx'' = -fir)df + ^ + r^dT.1, (8) 

where /(r) = r^/L^ + k - 2il//r, and rfS? = d9^ + 
siii^ {9)dip'^ with 6* G [0,7r], e [0, 27r] for the spheri- 
cally symmetric (k = 1), whereas dTjQ = d9^ + dip^ with 
0, S] — 00, oo[ for the planar case (k = 0) respectively. 
The plane symmetric black hole can also be seen as a par- 
ticular case of the spherically symmetric one, in the limit 
that the horizon is much larger than the AdS radius. 
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As previously discussed, spherically or planc- 
symraetric solutions of GR are solutions of DCS gravity. 
More precisely, for these highly symmetric spacetimes, 
*RR = = Cab- It follows that the SAdS spacetime is a 
solution to DCS gravity. This is also the case for black 
branes (sometimes called planar black holes) [13, 12^ . 
In fact, the theory reduces to Einstein-Klein Gordon 
theory with a (negative) cosmological constant and it 
is well known that such a theory does not admit scalar 
hairy black holes [l^ . As a consequence, the scalar field 
reduces to zero and we are left with pure GR with a 
cosmological constant. The situation can be different if a 
scalar self potential V{(f>) is introduced in the action [l8| . 
However in our case, even if the background solutions 
are the same as in GR, the dynamics of perturbations of 
such spacetimes can be modified due to the CS term. 



III. GRAVITATIONAL PERTURBATIONS OF 
SADS BLACK HOLES/BRANES IN DCS 
GRAVITY 



A. The equations 

Gravitational perturbations of Schwarzschild black 
holes were considered in Refs. [l^ [IB] in the context of 
DCS gravity in asymptotically flat spacetime. There, 
the authors showed that the metric perturbations fall in 
two categories: (i) a polar (even-parity) or Zerilli sec- 
tor, which is unaffected by the CS coupling, and whose 
equations of motion are identical to GR; and (ii) the axial 
(odd-parity) or Regge- Wheeler, which is modified in DCS 
gravity by the coupling to the dynamical scalar field. The 
procedure carries over to the case of asymptotically AdS 
spacetimes. In the following we focus only on the axial 
sector, since the polar sector is identical to the GR case 
(see e.g. Ref. [l3|)- We parametrize the metric perturba- 
tions around the SAdS background given by adopting 
the Regge- Wheeler gauge as in Ref. p^ . 

ds"^ = QabcLx^dx" + 2ho {r)Tdt + 2/ii (r)Tdr , (9) 

where 



T 



M0) 



d^YCje + F,{e)deY{;J^ , (10) 



and F^{e) = sin(6l) for k = 1, F^{e) 1 for k = 0, 
are the scalar harmonics on the sphere (plane) for k = 1 
(k = 0). Due to the background symmetry, we can fix 
m = without loss of generality. 

The perturbed gravitational equations then reduce to 



a set of three coupled equations: 

2rhi iuho , 2Mhi 



El 



L2 



/ 



0, 



1 



P ho{2f + {q-K)iq + 2K)) , 
E2 = ^fh„ + -lujfh. 



iufhi QMafa iMafa' 



r 



0, 



(w/iq iLoha hi (^r^Lj'^ — rf{q — K){q + 2k)') 



2/ rf 



2r3/ 



0, 



(11) 



where a — (j){r)r. Here and in the following, when k = 1, 
q is an integer [q ^ I = 2,3,...), whereas, for k = 0, 

q is real and we write q = ^Jqx + terms of the 

components {qx^ qy) of the spatial momentum on the 
plane, the scalar harmonics on the plane simply being 
^ eyi'p[i{qx9 + '7j,V')]- The three equations above are not 
independent and obey the following relation 



2ir^fE3 



2r^E2 i(q - K){q + 2k)Ei 



= 0. 



We supplement the equations by the derivative of Ei and 
solve the resulting set of four equations for hQ,h'Q, /iq and 
hi. Defining the Regge- Wheeler master function Q{r) as 



, rQ 
hi^ —p 



(12) 



the equations for Q and a reduce to 

dr^ 
f_ 
drl 

where is the tortoise coordinate defined as dr/dr^ 
fir) and 

6iMfa 



■^Q + {oj^ + VRw)Q + TRwa^O, (13) 
a + {uj^ + Vs)a + TsQ^O, (14) 



'q{q + k) 



1 



Trw 



2M 



2 

l2 



6m/M 

Ts (q + 2K)(q + K)q{q - k) g — 



(15) 



These equations form the basis of the present analysis 
and they reduce to the correct asymptotically fiat case 
[l^[l^ when L — > oo. Close to the horizon (r* — > — oo), 
the asymptotic solutions to this system are 



out iLor^ 



in —iLur.. 



(16) 

Close to spatial infinity (r* — )■ 0) on the other hand. 



wc have 



Ql 



oo 



irr 

oo 5 



(17) 
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Regular scalar perturbations should have a^" = 0, cor- 
responding to Dirichlet boundary conditions at infin- 
ity. The case for gravitational perturbations is less clear 
and the requirement of regularity is not unique. How- 
ever, most calculations in the literature assume Dirichlet 
boundary conditions also for the gravitational perturba- 
tions and this is the choice wc adopt here (see Ref. [14 1 
and references therein for a more detailed discussion 
on this problem). Therefore, we impose the following 
boundary conditions: 



Qout „out /niiT ^ii'i' n 



(18) 



In relation to the gauge/gravity correspondence, it has 
been shown that such boundary conditions allow one to 
recover the shear mode, but not the sound mode [l^, 
\2T\ . However, the sound mode is retrieved from the polar 
sector which is unaffected by the DCS coupling, and we 
therefore consider once and for all Dirichlet conditions at 
spatial infinity. 



B. Frobenius expansion 



In asymptotically flat spacetime, solving the coupled 
system analogous to Eqs. p^ - (fTi|) above is not an easy 
task [l^. However, in an asymptotically AdS spacetime, 
due to the boundary conditions (fT8|) . the coupled sys- 
tem p^ - ([T^ can be solved by extending the series ex- 
pansion method, used in Refs. [H, [2^ to study SAdS 
black holes in GR. 

I 



We start by re-writing / as 

^(,) ^ {r-r,){-L^+J^+rrn +rl) ^^^^ 

where r/j -I- r\/L^ = 2M. We factorize the horizon be- 
havior according to 



Q(r) = e-^"'--Q(r), a(r) 



•a(r), (20) 



and we introduce a new variable x = 1/r, rewriting the 
relevant equations p^ - p^ as 



(x-.,0.r(x)-;^ + A(.)A + !^ 
ax'' ax X — Xh 

cr a- 0, 

X~ Xh 

, X , X d'^ p , d vsix) 

(.x-..)..(.)^ + /.(.)^ + ^ 



Q{x) 



ts{x) 



Q(a;) =0, 



X- Xh 

where we have defined 
1 



x^f 



Xh ^ — , si{x) S2{x) = " •' , (21) 
rh X- Xh 

/i = /2 = 2x^f - x\f' - 2iiox^ , (22) 
vrw = {x - Xh)^^^^Y~ , = (x - x/i)— , (23) 



Vs Ts 
Vs = {x-Xh)—, ts = {x-Xh)y 

Expanding the unknown functions as 

oo / 1 \ ^ °^ 
Q{x)^^qkix j , a{x)^^ak 

k=0 ^ '^''^ k=0 

wc find the modified recursion relations 



(24) 



rh 



^ n — 1 

" k=0 
^ n— 1 

^ ~ p2 51 [('^(^ ~ l)s2,n-fc + fc/2,n-fc + Vs,n-k) CTk + ts^n-kQk] , 



k=0 



(25) 



where P,' = n{n — l)si,o + rifi.o and where we expanded 
the coefficients of the equation as 

oo 

X{x) = Y,Xk{x~Xh)\ (26) 

fc=0 

X generically denoting si , S2 , /i , /2 , , is , vrw , tRw ,o-,Q. 

There are two free parameters in the expansion: qq 
and (To J ^-nd all the coefficients may depend on w, so we 
have essentially (2 + 1) parameters. Either qo or ctq can 



be arbitrarily chosen due to the linearity of the equa- 
tions ([T ^ - p^ . Thus, we are left with two parameters, 
(To (or go) a-nd w, which can be found by simultaneously 
requiring Q{x — >■ 0) = and a{x — )- 0) = 0. The above 
procedure can be used to obtain the characteristic fre- 
quencies uj and the ratio ao/qQ. 
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IV. NUMERICAL RESULTS 

We solved the system of equations (|13p - P^ using the 
series expansion method described above. For numerical 
purposes, the series has to be truncated at some finite 
order A'^. As a first check on the numerical results, we 
verified that the corresponding solution a; at converges for 
large or moderately large N. A typical convergence plot 
is shown in Fig.|3]in the appendix, for two particular sets 
of data and for the spherically symmetric case. Usually 
the spherically symmetric case requires a larger trunca- 
tion order than that required for convergence in the pla- 
nar case. This is analogous to what happens in GR: the 
series method works well for large and intermediate black 
holes [13, [2^ , but it ceases to converge for small black 
holes. In the latter case, other methods (for example the 
resonance method [131) can be applied, but they do not 
extend naturally to coupled systems of equations as the 
one we deal with here. 

The truncation order needed to reach convergence 
strongly depends on the parameters. Roughly speaking, 
for a fixed accuracy the truncation order increases with 
a. For instance, our numerical experiments show that 
if K = 0, 30 — 50 to achieve an accuracy of a few 

tenths of a percent for any a/rf^ < 5. On the other 
hand, for spherical black holes with k = 1, sa 60 is 
required already for a/rf^ 1, this number increases for 
smaller spherical black holes. In what follows, we system- 
atically check consistency of the last points by increasing 
the truncation order. 

Note that, even if the limit where L — >■ 00 leads to 
the same system of equations than Ref. [l6j , the problem 
is not equivalent, due to the boundary conditions. As 
a consequence, it is not possible to recover results for 
the flat case. When a = 0, we recover the well-known 
modes of pure GR with a cosmological constant [3 [HI ■ 
Finally, in the hydrodynamic limit at small coupling, we 
found an analytic expression of the hydrodynamic mode, 
as described below. Our numerical results perfectly agree 
with the analytical result. 



A. QNM spectrum of AdS Schwarzschild black 
branes in DCS gravity 

Using the method described above, we computed the 
QNMs of a Schwarzschild black brane for different values 
of a and q. Our results arc summarized in Figs. [T] and 
[5) For axial gravitational perturbations and scalar per- 
turbations (these modes are coupled in DCS gravity) of 
Schwarzschild-AdS black holes, the fundamental modes 
are purely imaginary. In the left panel of Fig. [T]we show 
their dependence on a and on q. As expected, the CS 
corrections strongly depends on q. Indeed, modes with 
q = are plane-symmetric and the parity- violating CS 
contribution is vanishing in this case. In the right panel of 
Fig. [T] we also show the ratio (Jo/qo, which is purely imag- 
inary in this case. The ratio is vanishing when a = 0, 



signaling that this mode belongs to the gravitational sec- 
tor. In the large a limit the ratio asymptotically van- 
ishes, while it is finite for intermediate values of a. As 
discussed in Ref. [l^, the coupled system of equations 
can be understood in terms of a forced oscillator: the 
ratio cra/qo is related to the relative amplitude between 
the gravitational and the scalar waveform. 

In Fig. [5] we show the first overtones for gravitational 
and scalar QNMs. Here and in the following, the labels 
"scalar" and "gravitational" refer to the sector that is ex- 
cited in the GR {a — >■ 0) limit, but in DCS gravity both 
modes are proper oscillations of both gravitational and 
scalar perturbations. When a = 0, these modes are very 
similar, but they behave quite differently as functions of 
a. Indeed, the scalar mode is approximately constant, 
while the real part of the gravitational mode grows lin- 
early in the large a limit, while (although not shown in 
the plot) its imaginary part asymptotically vanishes. In 
Fig. [2] we show results for q = 3, but other values of q 
give the same qualitative results. Notice that in the large 
a limit the ratio Re[a;]/Im[a;] ^ 1 and the series method 
converges very slowly. Last points in Fig. [5] have been 
computed using iV ^ 80 as truncation order for the series. 
This prevents a tracking of the modes for larger values of 
a. Nevertheless, for some additional points, we increased 
the truncation order and we explicitly checked that the 
gravitational mode vanishes asymptotically. In particu- 
lar, in the region of the parameter space we spanned, i.e. 
up to a/rf^ < 1, none of the modes cross the real axis, 
i.e. we did not find any unstable mode also in the large 
a and large q limit. 



1. Hydrodynamic limit and AdS/CFT correspondence 

Having a large damping time, the purely imaginary 
mode dominates the thermalization timescale and the 
hydrodynamic limit of the 2 -f- 1 dual field theory [20| . 
Indeed, in the small-g limit, the location of this mode 
can be related to some hydrodynamic quantities (see for 
example Ref. [32|). 

For any c?— dimensional conformal field theory, the dis- 
persion relations for the shear mode and for the sound 
mode in the hydrodynamic limit read (ssj 

Wshcar = -iD^q^ - iD^iT^ + Tg)?* + O(g^) , (27) 

d-2 

Wsound ^ Vsq- i-, TD„q + 

a — 1 

d-2 d-2 \ ^ 

(28) 

where = 77/(6 + p) and 77, e, p and Vs are the shear 
viscosity, energy density, pressure and speed of sound, 
respectively, as computed from first order hydrodynam- 
ics. The dispersion relations above also get contribu- 
tions from higher order hydrodynamics: the relaxation 
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FIG. 1. Hydrodynamic mode for Schwarzschild black branes as a function of the DCS coupling ct, for different values of the 
momentum q. 




FIG. 2. First overtones for scalar and gravitational modes as functions of a for g = 3. The scalar mode is roughly independent 
from a. Different values of q give the same qualitative result. 



time is a second order quantity, whereas rs schemat- 
ically denotes all possible other contributions at order 
. In GR, these contributions come only from third- 
order hydrodynamics , which affect the shear mode at 
order j32| . However, in DCS gravity also first order hy- 
drodynamic quantities, like the Hall viscosity [isj . may 
affect the shear mode at the same order. 



On the other hand, the AdS/CFT duality relates the 
dispersion relations Eqs. and (HH) to the QNMs of 
SAdS black branes in the hydrodynamic limit. These 
frequencies can be computed analytically in the small q 
and small w limit. For example for four dimensional black 



branes we have [IJ, |26|, [2 

"2 L2 



^shcar 



^sound ^ i 



g - 1 

— I — ( 



(29) 



7 log 3 a/Stt 
9 2~ ^ ~2~ 



+ (30) 



where 7 = - 0.0281. Comparing Eqs. ^ 

and dSni) with Eqs. (H?]) and for d = 3, we can ex- 
tract the hydrodynamic quantities. More precisely, the 
coefficients of <^ in E qs. (1^^ and (|30p are related to the 
universal value of ?//s [ij, [21[ . The coefficient of q^ in the 
sound mode (|30p can be used to extract the relaxation 
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time in Eq. ([28]) . Finally, the coefficient of in the 
sound mode ([^ can be used to extract T3. 

For small CS coupling, we can extend the approxi- 
mate analytical methods used in Refs. [ssMsBl to our cou- 
pled system. First, we factorize the wavefunctions as in 
Eqs. (|20p . Then, anticipating the dispersion relation for 
the shear mode, we write 



(31) 



which holds in the hydrodynamic limit, g <C 1, and for 
small coupling. Finally, we solve the coupled system (|13|)- 
(|14p order by order imposing a series expansion of the 
following form: 

Q(r) = Qo(r) + QsW^^ + Q^{r)q^ + Q^{r)q^o? ,(32) 
a-(r) = aa(r)q^a. (33) 

It turns out that the ansatz above is consistent with 
the equations of motion and other possible terms (like 
those proportional to ~ aq^) vanish. The constants loi 
are found order by order by imposing Dirichlet boundary 
conditions on the solutions Qi (r) and Oi (r) [U • Remark- 
ably, in this limit the coupled equations can be solved 
analytically. At order aq"*, the scalar equation for a a is 
sourced by Qo c« 1/r. This can be solved analytically 
and its solution enters in the gravitational equation for 
Qa as a source term. The forms of Qiir) and ai(r) are 
cumbersome and not particularly interesting. Neverthe- 
less, the analytical expressions for Wi can be obtained in 
simple form. The first terms, loi and 0-14, are not affected 
by a and they are simply given in Eq. (|29|) . while the 
novel term oJq. reads 

= X « = — (201 - 20%/37r - 60 log 3) i - 0.123072i , 
640 V / 

(34) 

which is perfectly consistent with a fit of our numerical 
results shown in Fig. [TJ Therefore, the shear mode of a 
Schwarzschild black brane in DCS gravity, in the hydro- 
dynamic limit and for small CS coupling, reads 



3?'/i 



L2 



q^^O{q'>). (35) 



On the other hand, the sound mode belongs to the polar 
sector and it is not affected by the CS coupling. Our 
results imply that = (47r)~-^ also in DCS gravity. 
Moreover, since the sound sector is the same as in GR, 
also second-order hydrodynamic quantities like the relax- 
ation time (see e.g. Ref. [l^) are unaffected. The first 
non- vanishing correction appears in the q^ coefficient of 
the shear mode, which includes the Hall viscosity and 
other quantities related to second and third order hy- 
drodynamics. Comparing Eq. psp with the dispersion 
relations Eqs. (|27p and (pS)) , we obtain 



277 
4^ 



1 - 



X 
7L4' 



where we have set = 47rT/3, e + p 
relaxation time, 
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(37) 



(36) 



has the same value as in pure GR [3j]. 

Notice that, although in our case the background scalar 
field is vanishing, in our perturbed solution (j33p the 
scalar field ~' aq^ . This implies that the first non- 
vanishing correction to the Hall viscosity, r]A ^ cP'q^ (cf. 
Eq. (40) in Ref. [3), i-e. at the same order as T3. 

We remark that the universality of 77/s has been proved 
for theories without higher order derivatives. Indeed, in 
gravity theories with fourth or higher order derivatives 
7]/s is not universal |37l - [39} and the dual field theory has 
a different hydrodynamics, even at first order. Moreover, 
corrections to first order hydrodynamics can also come 
from Gauss-Bonnet terms in the gravitational action [iOl ■ 
Thus, it is remarkable that the CS coupling, although 
introducing higher curvature corrections, leaves 77/s and 
unchanged, while affecting only quantities related to 
the parity violation, entering the shear mode at order q^ 
in the hydrodynamic limit. 



V. CONCLUSION 

In this paper, we studied gravitational and scalar 
perturbations of planar and spherically symmetric 
Schwarzschild-AdS black holes embedded in dynamical 
Chern-Simons gravity. By generalizing a well-known se- 
ries expansion technique to compute characteristic fre- 
quencies in asymptotic AdS spacetimes to systems of cou- 
pled equations, we were able to study stability and oscil- 
lation modes of black hole and branes in this theory. We 
found that these solutions are stable also in DCS gravity, 
even in the large coupling and large momentum limit. In 
the opposite, hydrodynamic regime (small coupling and 
small momemntum), we explicitly computed the modifi- 
cation arising from the CS term to first relevant order in 
the couplings. By applying the gauge/gravity correspon- 
dence, we found that the entropy to viscosity ratio has 
the same "universal" value as in GR. Also the relaxation 
time is unaltered, while the CS coupling affects the shear 
mode at order g^, i.e. it affects the Hall viscosity and 
other second order and third order quantities related to 
parity violation. This constitutes one of the main results 
of this paper. It would be highly desirable to extend the 
analysis done in Ref. fisj in order to compute higher or- 
der hydrodynamic quantities. In such a way one could 
precisely relate the modifications (p6)) in the shear mode 
at order q^ to explicit hydrodynamic coefficients. 

One possible extension of this work is to study how 
the CS coupling affects the gravitational and scalar per- 
turbations of charged AdS black holes and the linear re- 
sponse of the holographic theory. One possible applica- 
tion in this context is, for example, including parity vi- 
olation terms in holographic superconductors [4l|. Also 
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in this case scalar and gravitational perturbations will 
be coupled through the CS coupling and this may affect 
the superconducting phase and the electrical conductiv- 
ity which, in the context of AdS/CFT correspondence, 
is computed from the Maxwell perturbations, coupled to 
gravitational perturbations. Again, the generalized series 
expansion can be used in this case. 

Another extension would be considering higher dimen- 
sions. For instance, in five dimensions, the CS term for 
gravity reads 

j (fytabcdeR'"'' A R^'' A e^d^x, (38) 

where i?"'' is the curvature form and e° is the fiinf-bein. 
In principle, this term should not be topological anymore, 
in the same way that the Gauss-Bonnet term is not a 
topological term in higher dimensions. It is naturally 
expected that spherically symmetric solutions to GR are 
still embeddable in DCS gravity. In this case, the dual 
theory is better known and the scalar operator in the 
dual theory can be explicitly identified. 



VI. ACKNOWLEDGEMENT 

This work was supported by the DyBHo- 
256667 ERG Starting and by FGT - Portu- 
gal through projects PTDG/FIS/098025/2008, 
PTDG/FIS/098032/2008, PTDG/GTE- 



AST/098034/2008 and GERN/FP/116341/2010. 
Computations were performed on the TeraGrid clusters 
TAGG Ranger and NIGS Kraken, the Milipeia cluster in 
Coimbra, Magerit in Madrid and LRZ in Munich. 

Appendix A: Spherical Black Hole 

For completeness, in this appendix we report our re- 
sults for the spherically symmetric black hole (k = 1). 
We investigated how the first overtone mode is changed 
by the CS coupling constant for small black holes, namely 
rh/L = 1 and q = 2. Fig. [3] shows how the results con- 
verge as function of the truncation order A'' as a increases, 
demonstrating a need for larger truncation numbers as 
the coupling increases. As we discussed for the planar 
case, our results exhibit a mixing between the the gravi- 
tationally sector to the scalar sector. In fact, both modes 
exist for a = where the equations are decoupled; in 
this sense, the modes are either purely gravitational, ei- 
ther purely scalar. Once a is set on, the modes start to 
mix, but in regimes where a^/qo ^ 1 (co/go ^ l)j the 
mode is dominated by the gravitational (scalar) sector. 
In this sense, the mode originating from the scalar sector 
in the a = theory stays dominated by the scalar sec- 
tor, whereas the mode originating from the gravitational 
sector is rapidly mixed with the scalar mode in the sense 
that the ratio cro/qo becomes of order 1. These consider- 
ations can be appreciated in Fig. |4l where the real and 
imaginary parts of the modes are shown, together with 
the ratio do/go- 
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